0^ 



(N 



> 



>$ 



Differences in the Cooling Behavior of Strange Quark Matter Stars and 

Neutron Stars 

Christoph Schaab^, Bernd Hermann, Fridolin Weber and Manfred K. Weigel 
Institut fiir theoretische Physik, Ludwig-Maximilians Universitat Miinchen, Theresienstr. 37, 

D-80333 Miinchen, Germany 

ABSTRACT 



^i \ The general statement that hypothetical strange (quark matter) stars cool more 

pLn I rapidly than neutron stars is investigated in greater detail. It is found that the direct 

Urea process could be forbidden not only in neutron stars but also in strange stars. In 
this case, strange stars are slowly cooling, and their surface temperatures are more or 
less indistinguishable from those of slowly cooling neutron stars. Furthermore the case 
of enhanced cooling is reinvestigated. It shows that strange stars cool significantly more 



00 . rapidly than neutron stars within the first ~ 30 years after birth. This feature could 

^ , become particularly interesting if continued observation of SN 1987A would reveal the 

O ' temperature of the possibly existing pulsar at its center. 

O,. Subject headings: stars: evolution - stars: neutron 

o 

52 ' 1. Introduction 



The theoretical possibility that strange quark matter - made up of roughly equal numbers of 

up, down and strange quarks - may be more stable than atomic nuclei (specifically iron, which is 

C^ I the most stable atomic nucleus) constitutes one of the most startling predictions of modern physics 

(see Bodmer 1971, Witten 1984, [Terazawa 1989| ), which, if true, would have implications of greatest 



importance for laboratory physics, cosmology, the early universe, its evolution to the present day, 
and massive astrophysical objects (cf. [Madsen fc Haensel 199l| ). Unfortunately it seems unlikely 
that lattice QCD calculations will be accurate enough in the foreseeable future to give a definitive 
prediction on the absolute stability of strange matter, so that one is presently left with experiments 
and astrophysical studies (cf. plendenning fc Weber 1992 , Glendenning et al. 1995 ) to either 



confirm or reject the absolute stability of strange matter. This letter, dealing with the second 
item, compares the cooling behavior of neutron stars with the one of their hypothetical strange 
counterpars - strange stars (cf. Witten 1984, Haensel et al. 19*8^ , Alcock et al. 1986| , plendenning 



1990). The theoretical predictions are compared with the body of observed data taken by ROSAT 
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and ASCA. There have been investigations on this topic prior to this one (e.g., see Alcock & Ohnto 



1988, Pizzochero 1991, Page 1992, Schaab et al. 1996b| ). These, however, did not incorporate the 



so-called standard cooling scenario that turns out to be possible not only in neutron star matter 
but in strange quark matter too, altering some of the conclusions made in the earlier investigations 
significantly. 



2. Description of strange matter 



We use the MIT bag model including 0(ac)-corrections (cf. phados et al. 1974 , Farhi &; 



Jaffe 1984) to model the properties of absolutely stable strange matter. Its equation of state and 



quark-lepton composition, which is governed by the conditions of chemical equilibrium and electric 
charge neutrality, is derived for that range of model parameters, that is, bag constant B^'^, the 
strange quark mass rUg, and strong coupling constant Oc, for which strange matter is absolutely 
stable (energy per baryon E/A less than the one in iron, E/A = 930 MeV). In the limiting case 
of vanishing quark mass, the electrons are not necesaary to maintain charge neutrality. In the 
realistic case of finite strange quark masse nis, the electrons can nevertheless vanish above some 



density which depends on ac- It was pointed out by Duncan et al. (1983) (see also Alcock et al. 



1986, Pethick 1992) that the neutrino emissivity of strange matter depends strongly on its electron 
fraction, Y^. For that reason we introduce two different, complementary parameter sets denoted 
SM-1 and SM-2 (see Tab. ffl), which correspond to strange matter that contains a relatively high 
electron fraction (SM-1), and 1^ = (SM-2) for the density range being of interest here. 



3. Neutrino emissivity 

The quark direct Urea processes 



and 



u + e + fe (1) 



u + e + i/e, (2) 



as well as their inverse ones are only possible if the fermi momenta of quarks and electrons (pp, 
i =u,d,s;e~) fulfill the so-called triangle inequality (e.g., pp < Pp+Pp for process (|l|)). This relation 
is the analogue to the triangle inequality established for nucleons and electrons in the nuclear matter 



case (direct Urea process, Boguta 1981 , Lattimer et al. 1991 ). 



If the electron fermi momentum is too small (i.e., Y^ is too little), then the triangle inequality 
for the above processes (|l]) and (|2|) cannot be fulfilled, and a bystander quark is needed to ensure 
energy and momentum conservation in the scattering process. The latter process is known as the 
quark modified Urea process, whose emissivity is considerably smaller than the emissivity of the 
direct Urea process. If the electron fraction vanishes entirely, as is the ease for SM-2, both the 



quark direct and the quark modified Urea processes become unimportant. The neutrino emission 
is then dominated by bremsstrahlung processes only, 

Q1 + Q2 — >Qi + Q2 + i' + '^, (3) 

where Qi, Q2 denote any pair of quark flavors. For the emissivities associated with the quark direct 
Urea, quark modified Urea, and quark bremsstrahlung processes, we refer to Refs. Price (1980), 
Iwamoto (1982) and Duncan et al. (1983). 



It has been suggested (see Bailin fc Love 197£ , Bailin fc Love 198^ ) that the quarks eventually 



may form Cooper pairs. This would suppress, as in the nuclear matter case, the neutrino emissivities 
by an exponential factor of exp(— A/A:bT), where A is the gap energy, k-Q Boltzmann's constant, 
and T the temperature. Unfortunately, there exists up to now neither a precise experimentally 
nor theoretically determined value of the gap energy. So to provide a feeling for the influence 
of a possibly superfluid behavior of the quarks in strange matter too, we choose A = 0.4 MeV 
as estimated in the work of Bailin & Love (1979). (Such a A value is not too different from the 
nuclear-matter case, where the proton ^Sq gap, for instance, amounts ~ 0.2—1.0 MeV (cf. [Wambach 



et al. 1991, Elgaroy et al. 1996), depending on the nucleon-nucleon interaction and the microscopic 



model.) The outcome of our superfluid strange matter calculations will be labeled SM-1 and 
SM-2^^ 



4. Observed data 

Among the X-ray observations of the 14 sources which were identified as pulsars, the ROSAT 
and ASCA observations of PSRs 0833-45 (Vela), 0656+14, 0630+18 (Geminga) and 1055-52 (see 
Tab. ^ achieved a sufficiently high photon flux such that the effective surface temperatures of 
these pulsars could be extracted by two- or three-component spectral fits (cf. pgelman 1995 ). The 



obtained effective surface temperatures, shown in Figs. || and ^, depend crucially on whether a 
hydrogen atmosphere is used or not. Since the photon flux measured solely in the X-ray energy 
band does not allow one to determine what kind of atmosphere one should use, we consider both 
the blackbody model and the hydrogen-atmosphere model, drawn in in Figs. |l| and |2| as solid and 
dashed error bars, respectively. The kind of atmosphere possessed by a specific pulsar could be 



determined by considering multiwavelength observations (see Pavlov et al. 1996). All error bars 
represent the 3(T error range due to the small photon fluxes. The pulsars' ages are determined by 
their spin-down times assuming a canonical value of 3 for the braking index. In reality the braking 
index may be quite different from 3. Its variation between 2 and 4, for instance, would change the 
age of Geminga as indicated by the horizontal error bar shown at the bottom of Figs. || and |2[ 



5. Results and Discussion 

The thermal evolution of strange stars and neutron stars was simulated using the evolutionary 
numerical code described in Schaab et al. (1996b) (see also Tsuruta 196(: , Richardson et al. 1982| , 



Van Riper 1991 , Page 1995 , Schaab et al. 1996a| ). The neutron star models are based on a 



broad collection of EOSs which comprises relativistic, fieldtheoretical equations of state as well 



as non-relativistic, Schroedinger-based ones (see Schaab et al. 1996b for details). As a specific 



feature of the relativistic models, they account for all baryon states that become populated in 
dense neutron star matter up to the highest densities reached in the cores of the heaviest neutron 
stars constructed from this collection of equations of state. Neutron stars are known to loose 
energy either via standard cooling or enhanced cooling. Both may be delayed by superfluidity. 
Consequently all four options have been taken into account here. These are labeled in Figs. |l| and 
as NS-1 (enhanced cooling) and NS-2 (standard cooling) for normal neutron star matter, and 
NS-1 and NS-2 (delayed cooling) for superfluid neutron star matter. The parameters of NS-l*^ 
and NS-2 are listed in Tab. 4 of Ref. Schaab et al. (1996b). In analogy to this, the corresponding 
strange-star cooling curves are SM-1 (enhanced cooling) and SM-2 (standard cooling) for normal 
strange quark matter, and SM-l*^^ and SM-2'^^ (delayed cooling) for superfluid quark matter. 

All calculations are performed for a star mass of M = IAMq, about which the observed 
pulsar masses tend to scatter. The band-like structure of the cooling curves is supposed to reflect 
the uncertainties inherent in the equation of state of neutron-star and strange-star matter. These 
have their origin, in the case of neutron stars (bands filled with dots), in the different many-body 
techniques used to solve the nuclear many-body problem, and the star's baryon-lepton composition. 
In the latter case, strange-star matter, the solid bands refer to different bag values, B^'^, which 
vary from 137 to 148 MeV for SM-1, and from 133 to 146 MeV for SM-2. All values correspond to 
absolutely stable strange matter. One might suspect that the large gap between the cooling tracks 
of the SM-1 and SM-2 models in Fig. |^ can be bridged steadily by varying the strong coupling 
constant ac in the range 0.1-0.15. However it turns out that the gap can be filled only for ac-values 
within an extremely small range. This is caused by the sensitive functional relationship between Oq 
and the neutrino luminosity L^, which is rather steep around that Oc-value for which the electrons 
vanish from the quark core of the star. All other values of Oc give cooling tracks which are close to 
the upper or lower bands, respectively. This behavior might be compared with the case of neutron 
stars, where the neutrino luminosity depends sensitively on the star's mass. 

One sees from Figs. ^ and ^ that, except for the first ~ 30 years of the lifetime of a newly born 
pulsar, both neutron stars and strange stars may show more or less the same cooling behavior, 
provided both types of stars are made up of either normal matter or superfluid matter. (We will 
come back to this issue below.) This is made possible by the fact that both standard cooling 
(NS-2) as well as enhanced cooling (NS-1) in neutron stars has its counterpart in strange stars 
too (SM-2 and SM-1, respectively). The point of time at which the surface temperature drop of 
a strange star occurs depends on the thickness of the nuclear crust that may envelope the strange 
matter core (cf. ^chaab et al. 1996b ). In the present calculation, strange stars possess the densest 



possible nuclear crust, which is about 0.2 km thick. Thinner crusts would lead to temperature 
drops at even earlier times, and thus an earlier onset of the photon cooling era. Figures |l] and ^ 
indicate that the cooling data of observed pulsars do not allow to decide about the true nature of 
the underlying collapsed star, that is, as to whether it is a strange star or a conventional neutron 
star. This could abruptly change with the observation of a very young pulsar shortly after its 
formation in a supernova explosion. In this case a prompt drop of the pulsar's temperature, say 



within the first 30 years after its formation, could offer a good signature of a strange star (see Alcock 



fc Olinto 1988 , Pizzochero 19911) . This feature, provided it withstands a rigorous future analysis 



of the microscopic properties of quark matter, could become particularly interesting if continued 
observation of SN 1987A would reveal the temperature of the possibly existing pulsar at its center. 

Finally, we add some comments about the possibility that only the neutron star is made up 
of superfluid matter but not the strange star. In this case one has to compare the models SM-1 
and SM-2 (see Fig. ||) with models NS-1 and NS-2 (see Fig. |^ yielding to an overall different 
cooling history of neutron stars and enhanced-cooling strange stars (SM-1). Therefore, the standard 
argument pointed out quite frequently in the literature that strange stars cool much more rapidly 
than neutron stars applies only to this special case. 



Table 1: Bag constant, B^'^, strange quark mass, mg, and the QCD coupling constant, Oc, for the 
two sets of parameters denoted SM-1 and SM-2. The energy per baryon, E/A, for two and three 
flavor quark matter is given too. 

Quantity SM-1 SM-2 



SV4 [MeV] 


140 


140 


TTis [MeV] 


150 


150 


Oc 


0.10 


0.15 


E/A in 2-flavor quark matter [MeV] 


959 


987 


E/A in 3-flavor quark matter [MeV] 


878 


892 



Table 2: Surface temperatures, T^, and spin-down ages, r, of observed pulsars. 



Pulsar 



T [yr] 



Model atmosphere 



T- [K] 



Reference 



0833-45 1.1 X 10^ blackbody 

(Vela) magnetic H-atmosphere 



1.3 X lO'' 
7.0tli X 10^ 



Ogelman et al. 1993 
Page et al. 199q 



0656±14 1.1 X 10^ 



blackbody 
magnetic H-atmosphere 



7.8^4:^ X 10^ 



Greiveldinger et al. 1996 
Anderson et al. 1993| 



0630±18 


3.2 X 10^ 


blackbody 
H-atmosphere 


5.2 ±3.0 X 10^ 
1.7 ±1.0 X 10^ 


Halpern & Ruderman 1993| 


(Geminga) 


|Meyer et al. 1994| 


1055-52 


5.4 X 10^ 


blackbody 


7.9lo-o X 10^ 


Greiveldinger et al. 1996 
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Fig. 1. — Cooling of non-superfluid strange star models SM-1 (lower solid band) and SM-2 (upper 
solid band), and neutron star models NS-1 (lower dotted band) and NS-2 (upper dotted band). 
The surface temperatures obtained with a blackbody- (magnetic hydrogen-) atmosphere are marked 
with solid (dashed) error bars (see Tab. |2|). The uncertainty in the pulsar's age is displayed by the 
error bar at the bottom. 

Fig. 2. — Cooling of superfluid strange star SM-l*^^ (lower solid band) and SM-2'^^ (upper solid 
band), and neutron star models NS-l*^^ (lower dotted band) and NS-2'^^ (upper dotted band). 
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